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ABSTRACT 

We present first measurements of the evolution of the scatter of the cosmic average early-type galaxy color- 
magnitude relation (CMR) from z = 1 to the present day, finding that it is consistent with models in which galaxies 
are constantly being added to the red sequence through truncation of star formation in blue cloud galaxies. We 
used a sample of over 700 red sequence, structurally-selected early-type galaxies (defined to have Sersic index 
> 2.5) with redshifts < z < 1 taken from the Extended Chandra Deep Field South (173 galaxies) and the Sloan 
Digital Sky Survey (550 galaxies), constructing rest-frame U — V colors accurate to < 0.04 mag. We find that the 
scatter of the CMR of cosmic average early-type galaxies is ~ 0. 1 mag in rest-frame U—V color at 0.05 < z < 0.75, 
and somewhat higher at z = 1 . We compared these observations with a model in which new red sequence galaxies 
are being constantly added at the rate required to match the observed number density evolution, and found that 
this model predicts the correct CMR scatter and its evolution. Furthermore, this model predicts approximately 
the correct number density of 'blue spheroids' — structurally early-type galaxies with blue colors — albeit with 
considerable model dependence. Thus, we conclude that both the evolution of the number density and colors 
of the early-type galaxy population paint a consistent picture in which the early-type galaxy population grows 
significantly between z = 1 and the present day through the quenching of star formation in blue cloud galaxies. 

Subject headings: galaxies: evolution — galaxies: general — galaxies: elliptical and lenticular — galaxies: 
stellar content — surveys 
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1. INTRODUCTION 

One of the best-known and most powerful scaling relations 
of the early-type (elliptical and lenticular) galaxy population is 
the systematic reddening of their colors with increasing lumi- 
nosity: the color-magnitude relation (CMR). The slope of the 
CMR is driven by a correlation between metallicity and mass 
( Faber & Jackson! 1 9761 i Kodama & Arimotol 19971 iTerlevich etail 
ll999tlTrager et alj200alGaUazzi et al.ll2006l) . while the scatter 
is determined by scatter in both age and metallicit y, where it is 



red sequence galaxy population (i.e., averaged over all envi- 
ronments) builds up in stellar mass by roughly a factor of two 
over the interval z = 1 to z = through the addition of new red 
sequence galaxies dChen et a l. 2003; Bell et alj|2004bl Icl matti 
l2006tlBrownet al. 2006at lFaber et al.l2007HScarlata et afT2 007, 
although some of the papers argue for a build-up in stellar 
mass only in galaxies with < 10 n Mo). These 'new' red se- 
quence galaxies are the result of truncation of st ar forma 



2000 argue that an anticorrelation between age and metallicity 



keeps the scatter of the early-type galaxy scaling relations rel- 
atively modest while allowing significant scatter in both age 
and m etallicity). The scatter in this correlation is relatively 
small (lBaumlll962HVisvanathan & Sandagel|l977l:lBower et all 
119921 ITerlevich et alj|200ll iMcIntosh et al.ll2005bl) . Because 
the color of stellar populations is strongly affected by their ages 
and metallicities (IWorth ev 1994), this correlation is a power- 
ful probe of the formation and evolution of the stellar popula- 
tions in early-type galaxies. The small intrinsic scatter found 
in so me clusters (as little as au-y = 0.04 ma g in Virgo and 
Coma; iBower et al. II 19921 ITerlevich et alj |200lh although other 
cluster s can have scatters approaching 0.1 mag; Mcintosh et all 
2005b) gave considerable momentum to the notion that early- 
type galaxies formed the bulk of their stars at early times and 
that their stellar populations have aged essentially passively to 
the present day. 

In the last few years, this position has been challenged by 
evidence from deep redshift surveys that the cosmic average 



generally thought that age is the dominant d river (IBower et al. _ 
1992 llTrager et alj2000llGallazzi et al.l2006l although lTrager et all 



tion i n some fraction of the blue cloud pop ulation jB ell et al. 
2007 ) through, e.g., galaxy-galaxy me rging dBell et al.l l2004b: 
Faber et all 120071: iHopkins et all 120071) or environmental pro- 



cesses such as strangu lation or ram-pressure stripping (e.g., 
iKodama & SmaifcOQll) . 

Such a scenario makes strong predictions about what the 
scatter o f the CMR and its evolution should be as a function of 
redshift (Ivan Dokkum & Franxll200ll) ; to first order the scatter 
is expected to be constant with redshift. The object of this paper 
is to quantitatively test this picture. We carefully measure the 
CMR evolution of galaxies selected to have concentrated light 
profiles and red colors — taken as a proxy for the early-type 
galaxy population — from z = 1 to the present day, and com- 
pare it to a toy model of a growing red sequence. We use the 
SDSS for the low-redshift CMR measurement, and a sample of 
galaxies with spectroscopic redshifts and accurate HST colors 
from the extended Chandra Deep Field South (CDFS hereafter) 
to probe the CMR out to z = 1 . As it was a priori unclear how 
large the CMR scatter should have been, we adopted a conser- 
vative approach that optimized rest-frame color accuracy at the 
expense of sample size, choosing galaxies with spectroscopic 
redshifts and within relatively narrow redshift slices (to mini- 
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FIG. 1. — The distribution of the spectroscopic redshift samples for the 
intermediate-redshift range. We used galaxies with redshifts around 0.55, 0.7 
and 1 .0 as indicated by the shaded areas, in order to minimize the contribution 
of ^-correction uncertainties to the CMR scatter. 



mize ^-correction uncertainties) and with color HST imaging 
(to minimize color measurement error). The data are described 
in §2. We describe the ^-corrections and their uncertainties for 
the intermediate redshift sample in §3. We present our mea- 
surements of the intercepts and scatter of the CMR for color 
and structurally selected samples as a function of redshift in 
§4. In §5, we compare the observed results to stellar popula- 
tion models of increasing complexity, finally comparing it to a 
model for the growth of the red sequence through the trunca- 
tion of star formation in blue cloud galaxies. We present our 
conclusions in §6. The casual reader may wish to skip to §5 
directly. In what follows, we assume £l,„.o = 0.3, £Ia ,o = 0-7, 
and Hp = 70 km s" 1 Mpc" 1 and a universally-applicable Kroupa 
(2001) stellar IMF for stellar population modeling. 

2. DATA 

We use data from two different sources, depending on red- 
shift. For galaxies with 0.5 < z < 1, we choose to determine 
accurate colors for galaxies with spectroscopic redshifts in the 
HST/Gems survey of the CDFS. This is crucial for minimiz- 
ing the color measurement error, placing the strongest possible 
constraints on the intrinsic scatter of the CMR. To complement 
this dataset at low redshift, we use a sample drawn from the 
SDSS at z = 0.05. 

2.1. Inte nnedia te redsh ift da ta 

For the intermediate redshift range, up to z = 1, we use data 
taken in the Extended Chandra Deep Field South (E-CDFS). A 
key ingredient is color HST imaging, taken from the Gems3 
and GOODS^ surveys. These data were used to estimate accu- 
rate galaxy colors within the half-light radius, and for selec- 
tion by galaxy structure. A second key ingredient, required to 
make precise fc-corrections, is accurate redshift measurements. 

1 Galaxy Evolution from Morphology and SEDs I Rix et aU2004f) 
2 Great Observatories Origins Deep Survey (Giavalisco et al. 2004) 



For this investigation, a sample of objects with reliable spec- 
troscopic redshifts was collected from a variety of sources. In- 
secure or uncertain spectroscopic redshifts were cross-checked 
with photometric redshifts from the COMBO - 1 ^3 survey. The 
spectroscopic selection criteria are described in the appendix 
and result in a final sample of 3440 galaxies with spectroscopic 
redshifts. We used 3030 galaxies in what follows; the remain- 
ing 410 objects could not be used either because they lacked 
data in one or more of the HST bands (211 galaxies), or because 
the object had a half-light radius smaller than 2 pixels (204 ob- 
jects and 5 objects satisfied both criteria). Figure [T] shows the 
redshift distribution of the 3030 useable objects. We choose the 
three best populated redshift bins of width Az = 0.1 at z ~ 0.55, 
0.7 and 1 .0 for our analysis. We made this selection for two rea- 
sons. Firstly, the bins have to be quite narrow to minimize the 
contribution of ^-correction (the wavelength range over which 
the interpolation must be done) uncertainties to the final error 
budget. Second, the best populated redshifts had to be used in 
order to obtain sufficiently populated CMDs for deriving CMR 
properties. In Figure Q] these redshift intervals are indicated by 
the shaded areas. 

For further analysis, we rejected known IR and X-ray sourcefl 
IR sources were identified by comparis on with MIPS 24/xm ob- 
servations of the CDFS from Spitzer dPapovich et alj2004 . We 
use the 24/im band owing to its high sensitivity to obscured 
star formation and AGN activity. The 80% completeness limit 
is 83/iJy, corresponding to approximate obscured SFR limits 
of (5, 10, 17) M^yr" 1 a t redshift s of 0.5 5 , 0.7 a nd 1.0 respec- 
tively dBell et al.l i2007b using a lKroupal (feoOl) IMF We fur- 
ther excluded galaxies detected in deep Chandra imaging of the 
CDFS. The coverage is non-uniform, with an exposure time of 
lMs in the central pointing, and 250ks per pointing in each of 4 
flanking fields. These depths are sufficient to detect moderate- 
luminosity AGNs (Lo.5-2keV = 10 41 — 10 42 ergs over the 
whole redshift range of interest dLehmer et alj|2005h . 

2.1.1. Colors and magnitudes from GEMS Images 

As the goal of this paper was to measure the scatter of the 
CMR as accurately as possible, we use high accuracy HST 
colors for the construction of the CMR. For all objects with 
spectroscopic redshifts we used Gems postage stamps to mea- 
sure accurate magnitudes in the F606W and F850LP passbands. 
We used the GALAPAGOS0 software package to cut postage 
stamps around the position of each object. Galaxy prop erties 
were adopted from a single-component Sersic model (Sersic 
1968) fit to th e 2-D galaxy lum inosity profile using the pack- 
age GALFIT dPeng et alj 120021) . We estimate colors within 
the observed half-light radius (not the intrinsic half-light ra- 
dius r e ,cALFir returned from fits to the light profile); such ob- 
served half-light radii are substantially larger than the intrinsic 
half-light radii for compact galaxies. We determined the ob- 
served half-light radius by performing aperture photometry in 
ellipses with the position angle and axis ratio given by GAL- 
FIT, out to 10 r f galfit on the F850LP-band image. The to- 
tal magnitude was defined as the aperture magnitude within 
10r e G4Lf7r- The observed half-light radius was then defined to 

3 Classifying Objec ts by Medium-Band Observations in 17 Filters 
IWolf et all20"oal2004l) 

4 Leaving in IR and X-ray sources gives similar results; we conservatively 
remove such sources in order to minimize the contribution of young stars and 
non-stellar light to the CMR scatter. 

5 Galaxy Analysis over Large Areas: Parameter Assessment by GALFITting 
Objects from SExtractor (M. Barden et al. in preparation) 
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FIG. 2.— The observed F606W-F850LP color within the half light radius 
plotted as a function of redshift. The red sequence and blue cloud are clearly 
visible. The cut between the two populations is schematically indicated by 
the dashed line (this cut is only for illustration purposes and not used in what 
follows for any analysis; the actual redshift and magnitude-dependent cuts are 
given in Sec. 4.1). The gray shaded areas show the three redshift ranges used 
in this paper. 



be the semi-major axis of the ellipse within which half of the to- 
tal light in F850LP was contained. In order to determine accu- 
rate F606W-F850LP colors, the F606W-image w as convolved 
with a difference PSF, determined from the PSFs dJahnke et al.l 
12004 in F606W and F850LP (the PSF correction was accurate 
to within 1 part in a million in terms of total flux on a pixel-by- 
pixel basis; see HauBler 2007 section 4.2.2 for details), and the 
flux in F606W within the F850LP half-light ellipse was mea- 
sured. Figure [2] shows the distribution of the measured colors 
against redshift. Uncertainties in the measured fluxes include 
contributions from Poisson uncertainty, read noise, and uncer- 
tainty from inaccuracies in the assumed sky level (this contri- 
bution is equal to the half-light area in pixels times the sky level 
uncertainty in counts per pixel). Typical values for uncertainties 
in the determined magnitudes are around 0.007 (slightly higher 
for high redshifts and lower for smaller redshifts). 

2.2. The low-redshift sample from the Sloan Digital Sky 

Survey 

In order to explore the scatter of the CMR at low redshift, we 
use a sample of galaxies draw n from the Sloan Digital Sk y Sur- 
vey (SDSS) Data Release 4 d Adelman-McCarthvl 120061) . We 
choose a sample of galaxies from the publicly-available New 
York University Va lue-Added Galaxy Catalog (NYU VAGC; 
Blanton_etal]L2005) in a very narrow redshift range 0.0495 < 
Z < 0.0505 (2053 galaxies). To maximize the accuracy of the 
color information, we adopt colors derived from 'Model' mag- 
nitudes — such magnitudes use the best-fit de Vaucouleurs 
or exponential models in the r-band as a kernel for measur- 
ing fluxes in ugi^ We fc-correct the observed-frame ugriz 
model colors to rest-frame U — V color using InterRest (the 

6 We also used aperture magnitudes with aperture radii 5 < r/arcsec < 10, 
finding similar or larger CMR scatter. 



FIG. 3. — An assessment of the ^-correction uncertainties through compar- 
ison of the observed F775W magnitude of those sample galaxies that overlap 
with the GOODS field and the F775W magnitude predicted using InterRest us- 
ing only the redshift, F606W and F850LP magnitudes as input. In this plot we 
show the difference between measured and calculated magnitude in F775W as 
a function of redshift. The diamond symbols indicate red objects with a Sersic 
index n > 2.5. The outlier at z ~ 0.1 is a composite merging system with both 
a dust-reddened and a blue component; such a composite system is too com- 
plicated to be described by the simple templates used here. The long dashed 
line shows the mean of these values and the short dashed lines show the RMS 
area around the mean value (both values are around 0.03). The random error 
serves as an approximation of the ir-correction uncertainty. 

same method we use for the intermediate redshift data; see Sec- 
tion[3]). Rest-frame magnit udes were calc ulated by ^-correcting 
Sersic model magnitudes (Blanton et aTl l2003f) from the NYU 
VAGC; such magnitudes are closer to the total magnitudes than 
Petrosian or Model magnitudes. The formal random error in the 
final U—V color is < 0.02 mag from photometric error, with an 
estimated calibration error of < 0.04 mag (assuming 0.01 mag 
calibration errors in g and r bands); random errors in My are 
^0.1 mag, dominated by systematic errors in how one defines 
sky levels and total magnitude. 

3. ^-CORRECTIONS 

To compare measurements in a redshift-independent manner, 
we A:-correct the observed frame measurements into restframe 
properties (throughout this paper we use My. rest and (U - V) rest 
although our conclusions do not depend on this choice). To do 
so we use the IDL implemented restframe interpolation code In- 
terRest by ENT (http://www.strw.leidenuniv.nl/^ent/InterRest; 
this work is based on an earlier version of the code). To de- 
rive the redshift-dependent transformation between observed 
and rest-frame colors (t he algorithm is de scribed in more de- 
tail in Appendix C of iRudnick et alj |2003). InterRest uses ob- 
served SEDs from a set of templa t e gala xies: four empirical 
model spectra from Colema n et al.l dl9 80). and one additional 
starburst template from lKinnev et al.l ([1996); this is a subset of 
the default template set which is used to avoid degeneracies in 
color space. 

These ^-corrections were tested in two ways. Firstly, one 
can use the ^-correction routine to predict a F775W magni- 
tude for galaxies in GOODS (where one has F606W, F775W 
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and F850LP). We show such a comparison in Fig. [3] where one 
can see that the ^-corrections give a scatter of ~ 0.03 and an 
equal amount of mean color offset. Secondly, we can com- 
pare the final rest-frame colors to the results of independent 
^-correction codes. Comparison to (lower accuracy) COMBO- 
17 colors shows a scatter of 0.12 magnitudes and an offset 
of 0.06 mag, and comparison to a stellar population model- 
derived ^-correction by EFB and BH shows measurement-to- 
measurement scatter of less than 0.02 mag but overall rest- 
frame color offsets of ^0.1 mag (this results primarily from 
a difference between the 4000A break structure of the stellar 
population models and the observed ones used by InterRest). 
We conclude that the ^-corrections of the HST-derived colors 
are accurate to a few hundredths of a magnitude (random error; 
this contributes to the scatter on the CMR; Taylor et al., 2009, 
submitted), with possible overall systematics of < 0.1 mag (af- 
fecting primarily the zero-point of the CMR). We assume the k- 
correction errors to scale with the wavelength range over which 
the interpolation is done. 

4. RESULTS 

In this section we present our analysis of the red sequence 
CMR scatter and its evolution. First, we explain the selection 
criteria applied for red sequence galaxies, then we discuss the 
actual measurement of the scatter of the red sequence at differ- 
ent redshifts 0.05 < z < 1.05. With our precise magnitudes and 
colors we have produced color-magnitude diagrams (CMDs) 
for four redshift intervals and plot them in Figure|4] 

4. 1 . Selection of Red Sequence Galaxies 

Our selection criteria take care to choose only non-star- 
forming, early-type galaxies representing the red sequence and 
to reject blue galaxies, galaxies with disk structure, and galax- 
ies with ongoing star formation and/or AGN activity. 

First, we applied a color cut to exclude galaxies with obvious 
signs of star formation from the sample. For this we made use 
of the bimodality of the color distribution which is visible in 
the CMDs at all redshifts (Figure|4|i. The color cut was chosen 
to fall into the gap between blue and red objects (see also his- 
tograms in Fig. HJ. We used a tilted cut with the same slope 
for all samples. As the mean colors of the populations evolve 
with redshift we take this into account by allowing the cut to 
evolve correspondingly. The cut applied here can be described 
as (U-V)rest > -0.085 -My^-0.65-0.5 -z. Small changes in 
the positions of the cuts do not influence the general results. In 
the diagrams the cut is made visible by different shades of gray. 
The light gray points are objects bluer then the color criterion 
whereas the dark gray and black points indicate object on the 
red side of the cut. 

A second selection criterion was applied to clean the sam- 
ple of structurally late-type galaxies. For this the Sersic in- 
dex of the galaxy is used (for the Gems galaxies, this is mea- 
sured using Ga lfit, while for the S PSS galaxies we adopt the 
Sersic fits from Blanton et al. I I2005I) . To be treated as part of 
the red sequence a Sersic index of n > 2.5 is required. This 
cut successfully weeds out edge-on spirals from the sample 
dMcIntosh et al. 2005a). These objects are shown as dark gray 
points. Regardless of color and Sersic index we reject galaxies 
with X-ray or 24/im detections (for the CDFS data only; the 
SDSS sample we use lacks such information) as their colors 
may be substantially affected by UV-bright young stellar pop- 
ulations and/or an accreting supermassive black hole (although 



our results remain unchanged if these systems are not excluded 
from our analysis). These objects are shown in the CMDs with 
asterisks (X-ray) and diamonds (IR). The objects which qual- 
ified as red sequence galaxies after these selection criteria are 
shown in black with error bars in the CMDs (number of RS 
galaxies in the redshift bins: 563/38/102/31). 

4.2. Fitting the Red Sequence 

The goal of this paper is to measure the intercept of the CMR 
at a fiducial magnitude and its scatter, and to show the evo- 
lution of both the CMR intercept and scatter as a function of 
redshift (see Figure 0J. The intercept is the color value of the 
relation measured at a specific magnitude (either the same mag- 
nitude value at all redshifts, or using a variable magnitude that 
attempts to account for the evolution of the population). The 
scatter is calculated using the IDL routine robust _sigma, which 
gives an outlier-resistant estimate of the dispersion of the color 
distribution using the median absolute deviation as an initial 
estimate and then weights points using Tukey's Biweight as an 
robust esti mator (this is the usual approach employ ed in this 
field; e.g.. iBower etaT]|l992b iMcIntosh et al.ll2005bl) . We fit 
a simple linear function to the distribution of colors and mag- 
nitudes of the non-star-forming early-type galaxies under two 
assumptions: 

1 . We fit the data points in each redshift bin separately. 

2. We choose to hold the slope fixed at all redshifts > 0.5 
to the value measured at z = 0.05 from the SDSS. 

The slope of the two methods are in some cases very differ- 
ent, but as the SDSS sample has many more objects and a 
larger dynamic range, this fit is much more reliable than the 
fits for the other bins. Further support for this assumed slope is 
given by the best populated Gems bin around z = 0.7, which 
shows nearly the same slope. We note that the determina- 
tion of the intercept and scatter are not sensitive to our choice 
of fitting method (Table HB The histogram of color offsets 
from the CMR with redshift-independent slope is shown in the 
inset panels of Fig. [4] for early-type galaxies (black) and for 
all galaxies (gray). As the observed scatter contains contri- 
butions from intrinsic scatter and measurement uncertainty we 
estimate the intrinsic color scatter by subtracting the random 
color uncertainties in quadrature. The color uncertainties for 
the low redshift galaxies were estimated by scaling the Model 
u — r color uncertainties by d(U — V) res , jd(u — r), and adding 
in quadrature a small empirically-determined 0.015 mag con- 
tribution accounting for real scatter in U — V at a given u-r. 

S(U-V) 2 = [S(u-r)-d(U-V)/d(u-r)] 2 + 0.015 2 . If not fur- 
ther specified we refer to the intrinsic scatter hereafter. 

The color intercepts (U-V) res , of the fitted CMR at a fixed 
Mv,rest = _ 22 are shown as a function of redshift in Figure [5] 
The solid errorbars denote the formal error in the mean inter- 
cept, calculated using bootstrapping, although recall that the 
actual, primarily systematic, uncertainty in the galaxy colors is 

furthermore the effect of excluding IR and X-ray detections from the red 
sequence is quite small and does not affect our results. The effect on the inter- 
cepts is around 1%. The values for intrinsic scatter would be a bit larger if IR 
and X-ray sources are included in the sample, changing for z = 0.55 from 0. 128 
to 0.138, for z = 0.7 from 0.126 to 0.133 and for z= 1.0 from 0.207 to 0.216. 

8 Recall that U - V is estimated from the ugriz. Model colors using InterRest. 
Our choice of scaling from u — r uncertainty is reasonable: u, g and r are the 
main determinats of U — V , where g and r are relatively well-measured and the 
u band uncertainty is factors of several larger. 
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FIG. 4. — U — V restframe color against My for four different redshift bins. The light gray points are objects below a tilted cut, dividing blue and red objects. 
This cut has a fixed slope and a redshift-dependent intercept. Depending on the Sersic index the red objects are divided in a dark gray (n < 2.5) and a black sample 
(n > 2.5; with errorbars for the three intermediate redshift bins); the n > 2.5 red galaxies are retained in our analysis. Independent of color and Sersic index IR 
(diamonds) and X-ray (asterisks) detections are rejected from the red sequence sample. The dashed lines are fits to the red sequence samples for the three redshift 
bins with z > 0.5. In the z = 0.05 panel the fit is shown as a solid line (the position of this fit is shown in the other panels as the dotted line). A line with the same 
slope is shifted to the mean color value of the three other panels (solid line). Using these 'fits' we calculated intercepts at My = —22. The small histograms show 
the color distribution measured as an offset from the relation. The complete distribution is shown in light gray, whereas the red sequence (the black points in the 
CMDs) is shown in black. The dashed vertical line shows the intercept of the fit (dashed line in the CMDs) at a magnitude value of —22. The solid vertical line 
shows the mean intercepts with the CMR scatter (RMS) marked by the dotted vertical lines. The dash-dotted lines show the CMR intercept at z. = 0.05. Overplotted 
is a Gaussian with mean and <r defined by the mean color and RMS of the red sequence sample. 



^0.1 mag. The dotted 'errorbars' show our measurement of the ments from othe r studies dBlantonll2006t iFranzetti et al.ll2007l 

CMR scatter, calculated as the resistant dispersion of the color ICool et al]|2006|) are shown in this plot together with our data, 

offsets of the early-type galaxies from the CMR. The results In the case of lBlantonl d2006l) and lCool et all (120061) . we have 

are tabulated in Table Q] For comparison, intercept measure- ^-corrected their values to rest-frame U — V color using the 
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Table 1 









Fit 


measured 


intrinsic 




Fixed slope 


measured 


intrinsic 






slope 


intercept 


scatter 


scatter 


slope 


intercept 


scatter 


scatter 


z = 


0.05 


-0.087 ± 0.006 


1.457 ± 0.005 


0.124 


0.117 


-0.087 


1.457 ±0.005 


0.124 


0.117 


z = 


0.55 


-0.015 ± 0.028 


1.206 ± 0.021 


0.126 


0.124 




1.248 ±0.022 


0.130 


0.129 


z = 


0.70 


-0.066 ± 0.014 


1.178 ±0.010 


0.127 


0.127 




1.184 ±0.011 


0.126 


0.126 


z = 


1.00 


0.010 ± 0.047 


1.097 ± 0.041 


0.210 


0.206 




1.103 ±0.036 


0.212 


0.208 



Note. — Color intercepts (at My.rest = —22) and scatter in (U — V) res , measured with two different methods. 'Fixed slope' means that the slope of the CMR is fixed 
to the z = 0.05 value for all redshifts (solid lines in Fig. [4), while in the 'Fit' columns the red sequence is fitted for each bin (dashed lines in Fig. [4)- 
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FIG. 5. — CMR intercepts in U -V as a function of redshift. The small solid 
errorbars denote the standard deviation of all intercept values calculated with 
the bootstrap method. The 'errorbars' with dotted lines indicate the resistant 
dispersion of the color offsets from the best-fit CMR (i.e., the scatter of the 
CMR). The other symbols show measurements of the CMR intercepts from 
other papers, as explained in the legend. The measurements from Blanton and 
Cool were transformed into U — V res , using kcorrect. The errorbars reflect the 
uncertainties introduced by this treatment. 



kcorrect software package, written by M. Blanton. We made 
sure that there is no significant difference (< 0.03 mag) between 
the rest-frame colors derived with kcorrect and InterRest. 

We find that the intercept of the color-magnitude relation 
evolves by ~ 0.3 mag in U — V color between z = 1 and z = 
0.05, in rough agreement with previous studies (e.g. | Bell et alj 
I2004at iFranzetti et alj|2007h lBlantonl l2006; Tay lor et alj |2008), 
although recall that systematic errors in the intercepts of the 
CMR are significant, ~ 0.1 mag in U — V rest-frame color. 
More importantly, we find a scatter of <Ju-v ~ 0.1 mag (Table 
[TJ at all z < 0.75; the scatter at z = 1 appears t o be somewhat 
larger (being consistent with lTavlor et aLlfeOOSt . These values 
are somewhat larger than the scatter measured for local galaxy 
clusters (0.04 < au-v < 0.1 , where the scatter app ears to vary 
from cluster to cluster; e.g. iMcIntosh et al. 2005b£J We note 

'The possible trend towards a bluer CMR intercept at very large cluster- 
centric radii, arguably an environmental effect, would also incr ease the scatter 
of the cosmic ave rage CMR relative to the cluster CMR (e.g.. Terlevi ch et all 
l200lt rPimbblet et al. 200%). 



that subtraction in quadrature of the individual galaxy measure- 
ment errors from the CMR scatter yields almost unchanged re- 
sults, even for the low redshift SDSS sample (where the scatter 
could decrease to ~ 0.09 mag). This shows that we measured 
a real scatter in the relation and not only a spread caused by 
measurement uncertainties. In particular, one should note that 
for the SDSS, spectral analysis of the drivers of the CMR scat- 
ter have demonstrated clear spectral differences between early- 
type galaxies of a given magnitude at the red side and blue side 
of the CMR dGallazzi et all 120061; ICool et alj |2006l), demon- 
strating that the bulk of the CMR scatter is intrinsic. 

5. INTERPRETATION 

In the last sections, we described our measurements of the 
CMR intercept and scatter for cosmic average early-type galax- 
ies at four redshifts between 1.0 and 0.05. We find a slowly- 
evolving mean color, and an almost non-evolving scatter of 
^0.1 mag in U-V rest-frame color. In this section, we build 
some intuition about possible interpretation of this result us- 
ing stellar population synthesis modeling. We show first the 
evolution of single bursts, or populations of galaxies whose star 
formation is truncated at a particular time, in §5. 1 and 5.2. Note 
that the purpose here is not to test 'monolithic collapse' models 
(those models are already ruled out by the observed build up of 
the red sequence population), rather it is to establish the basic 
model ingredients. In §5.3 and 5.4 we ask and answer a simple 
question (similar in spirit to lvan Dokkum & Franxl200lt) : if the 
red sequence forms through ongoing and continious truncation 
of star formation in previously star-forming systems, can one 
simultaneously reproduce the build-up of red sequence galax- 
ies and the scatter of the CMR? 

5.1. The evolution of single bursts 

To get our bearings, we focus first on the evolution of the 
CMR intercept using single bursts of star formation. In what 
follows, we use the stellar population model P egaseFI (Ver- 
sion 2.0, see iFioc & Rocca-V olmerangel l 19971 for the descrip- 
tion of an earlier version of the model) to predict the evolution 
of galaxy colors and absolute magnitudes as a function of red- 
shift. 

The results are shown in Figure [6] for three different forma- 
tion redshifts and for various metallicities. One can see that 
the evolution of the intercept of the CMR follows roughly the 
trends expected for the passive reddening and fading of ancient 
stellar populations. In the context of single bursts, the rate of 
change of the CMR intercept with redshift is related primarily 

10 Projet d'Etude des GAlaxies par Synthese Evolutive 
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FIG. 6. — The evolution in rest-frame color of single-burst stellar populations formed at three different redshifts, and with different metallicities (where Z = 0.02 is 
solar metallicity). The asterisks denote the measured intercepts from our paper; the gray asterisks are the same as in Fig. \5\ whereas the black asterisks show color 
intercepts measured at a redshift-dependent absolute magnitude (to compensate for fading of stellar populations as the galaxy population ages; one can think of this 
as measuring the intercept at approximately constant stellar mass). The other gray symbols show data from other papers, measured in a similar fashion to our gray 
data points (see the legends of Fig. [5]or[7}. 



to formation history, whereas the overall intercept is sensitive 
primarily to metallicity, with some sensitivity to age. 

To compensate for the passive fading of the early-type pop- 
ulation, we also measure the CMR intercept at a redshift- 
dependent absolute magnitude (black asterisks in Figure[6]). We 
choose to adopt a model with metallicity Z = 0.008 and a for- 
mation redshift of 2 to estimate the evolution of the luminosity 
of an early-type galaxy as a function of redshiff^l and use this 
to calculate the absolute magnitude values at which we measure 
the intercept of the CMR (these absolute magnitudes are —22, 
-21.5, -21.3, and -20.8 at z = 1.0, 0.75, 0.55 and 0.05 respec- 
tively). One can think of these intercepts as being measured at 
a given, redshift-independent, stellar mass. Because these in- 
tercepts reflect the evolution at a given stellar mass, we choose 
to focus on these values in what follows, as a more intuitive re- 
flection of the likely evolution of the color of a given galaxy (or 
a mass/metallicity bin in the evolving galaxy population). 

5.2. The evolution of galaxies with truncated star formation 

Motivated by the observational evidence for truncation of 
star formation in blue sequence galaxies, leading to the ongo- 
ing build-up of the red sequence galaxies at redshifts z < 1 (e.g. 
iBell et alJl2004bt iBorch et alll200a [BelTeraD[2007h . we study 
the evolution of the colors of galaxies whose star formation has 
been truncated at a variety of redshifts z < 1 . The purpose of 
this section is to get a feeling for the ingredients used for the 
modeling of the CMR evolution described in the next section. 
Figure [7] shows the rest-frame color evolution of a stellar pop- 
ulation which forms stars at a constant rate from z/ • = 2 to the 
truncation redshift z t runc- We plot the evolution of seven dif- 
ferent truncation histories which all have a constant metallic- 
ity of Z = 0.02. The rightmost track has a truncation redshift 
very close to the formation redshift of the model, namely at 
z = 1.95. The other lines show the evolution of galaxies with 
smaller truncation redshifts (from left to right) from 0.2 to 1 .2 

1 1 The choice between the formation redshifts in the range shown in the plot 
does not have a great influence on the results; to choose z/ = 3 instead would 
lead to in a slightly smaller evolution in magnitudes, but the resulting difference 
in intercepts is smaller than 0.02 for all redshifts. A different choice in metal- 
licity has also only very small effects. To choose Z = 0.02 (solar metallicity) 
would change the results even less than the change of formation redshift. 
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FIG. 7. — The redshift evolution of rest-frame color for different truncation 
times. For all lines, star formation starts at z = 2 and the different lines denote 
a variety of different truncation redshifts (from left to right: 0.2 to 1.2 in steps 
of 0.2). The tiny black dots on the lines indicate the point on the evolutionary 
track a galaxy has reached 0.5, 1 and 2 Gyr, respectively, after the truncation 
of star formation (as labeled on the second line from the right). The symbols 
are the same as in Figures[5]and|6] 



in steps of 0.2. A general feature of such truncation models 
is a period of relatively rapid evolution onto th e red sequence 
(timescales < 1 Gyr, as can be seen in Fig.Hand|Hintoil2006; 
ISchweizer & Seitzerl 19921) . and then relatively slow subsequent 
fading and reddening of the population. 

5.3. Expectations for the Scatter Measurement 

With these SF modeling ingredients in place, we are now in a 
position to ask what the evolution of the intercept and scatter of 
the CMR can tell us about the evolution of the early-type galaxy 
population. Obviously, we will be unable to address this issue 
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FIG. 8. — The integrated number of truncated (pre-red sequence) galaxies as 
a function of cosmic age (t ga / lnmc = 13.5 Gyr is the present day). Assuming 
that after truncation of star formation galaxies fade onto the red sequence in 1 
Gyr, one reproduces (by construction) the observed evolution in stellar mass 
density on the red sequence (inset panel). 
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FIG. 10. — Colors of individual galaxies calculated in the same way as for 
Figure|9]but assuming a scatter in metallicity. The lines show the mean value 
in color and the scatter of the color distribution. Recall t hat w e have used a 
redshift-dependent luminosity cut as described in Section [5~71 to measure the 
intercepts. 



completely, as we have seen that many factors influence the 
colors and magnitudes of galaxies: when star formation started, 
metallicity, when (if!) star formation ends, and whether or not 
there is any residual star formation. 

Instead, we limit ourselves to one well-defined question. Re- 
cent observations have measured a significant increase in the 
total s tellar mass in the red sequence galaxy population since 
z = 1 dChen et al.l l2003t iBell et alj l2004bt iBundv et all 120051; 
iBrown et al.ll2006bl; iBorch et al.ll2006t iFaber et alj|2007l) . This 
evolution manifests itself primarily in terms of a n increasing 
space density of red s equence galaxies at < L* dBorch et al.l 
2006l: lFaber et alJ20 07). and is most naturally interpreted as be- 
ing fed by the trunc ation or quenching of star formation in mas - 
sive blue galaxies dBell et alj|2004bl l2007t IFaber et al.ll2007l) . 
In such a scenario, one expects a significant scatter in the CMR, 
because of the constant flow of recently star-forming galaxies 
onto the red sequence. Here, we will predict the CMR scat- 
ter implied by a constantly growing red sequence, and compare 
this predicted scatt er with the observations. Su ch an exercise 
was carried out by Ivan D okkum & Franxl d200ll) for early-type 
galaxies in clusters; here, we extend their work to a cosmic av- 
erage environment. 

Our first ingredient is a toy model for the build up in the 
number of red sequence galaxies as a function of time. We 
assume i) that the growth in the total mass on the red sequence 
is driven entirely by adding blue galaxies that have had their 
SF recently truncated, and ii) red sequence galaxies are added 
at all stellar masses equal hPl We describe the evolution of 
the number of red sequence galaxies with a simple linear fit 
(in redshift) as shown in the small panel of Fig. [8] compared 
w ith the meas urements of the integrated stellar mass density 
of IBorch et alj d2006l) . We then determine the derivative of this 
relation to get the change in stellar mass density with time. This 

12 In fact the CMR scatter and mass build-up are both measured at ~ 
10 u Mq, reducing the importance of this assumption. 



is used to estimate the truncation rate by shifting the time axis 
by -1 Gyr - approximately the time ta ken to redden e nough to 
satisfy our red sequence cut (following Blanton 2006; see also 
Fig. [7]). In thi s way, we have a truncation rate history that is 
consistent with Borch et al. I d2006h and can be used to estimate 
the evolution of the CMR sc atter. To describe the truncation 
rate for epochs not covered bv lBorch et alj d2006l) . we choose to 
model the truncation rate as a constantly increasing function of 
time fromz/ ~ 5.5 until z trunc = 0.9. The influence of this choice 
on the results is negligible. The resulting number of truncated 
(pre-red sequence) galaxies - the integral of the truncation rate 
- as a function of cosmic epoch is also given in Fig. [8] 

The next step is to build a more realistic truncation model 
(c.f. Fig. |7): instead of choosing one particular truncation time 
as in §5.2, we instead draw truncation times from the truncation 
time distribution in Fig. [8] We choose 12 'observation' red- 
shifts between z = and z = 1.1; at each we draw 1000 galaxies 
from this truncatiorPl history shown in Fig.[8]resulting in sam- 
ples of galaxies with different truncation redshifts, but the same 
formation redshift z s tart = 3 for all galaxieQ Due to this variety 
of truncation redshifts the galaxies in the samples experienced 
periods of passive evolution of different length. This then gives 
a distribution of galaxy colors at each redshift of interest; for 
display purposes in Figs.l9land[T0lwe have added a small ran- 
dom offset in redshift. We assume no ongoing low-level star 
formation in red sequence galaxies; while potentially unrealis- 
tic, it allows us to estimate the expected scatter from truncation 

13 The results of this paper do not change significantly if a burst of star forma- 
tion occurs before truncation, as might be expected in a gas-rich galaxy merger. 

14 Note that this choice of Zstart differs from zr in section f5~Tl In contrast to 
this former Zf, which was used to simulate the passive evolution after an initial 
star burst, Zsian indicates the starting point of a longer period of star formation. 
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FIG. 9. — Colors of individual galaxies calculated using the truncation history shown in Fig. fusing one fixed metallicity for each plot. The colors were measured 
for certain redshift values and are shown here spread out in redshift to illustrate the density in color space. The metallicities are from left to right 0.008, 0.02 (solar) 
and 0.05. There is a bimodality at some redshifts, as a result of kinks in the evolutionary tracks at some values of age/metallicity. 
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Fit 




Fixed slope 




Simulation 




measured 


intrinsic 


measured 


intrinsic 






scatter 


scatter 


scatter 


scatter 




z = 0.05 


0.124 


0.117 


0.124 


0.117 


0.115 


z = 0.55 


0.126 


0.124 


0.130 


0.129 


0.126 


z = 0.70 


0.127 


0.127 


0.126 


0.126 


0.132 


z= 1.00 


0.210 


0.206 


0.212 


0.208 


0.167 



Note. — For better comparison the measurement values for the measured and intrinsic red sequence scatter are presented here together with the values for the 
simulated scatter. 
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FIG. 1 1 . — Metallicity distribution for the modeled galaxy sample in Figure 
1101 The mean value (ZogioZ/Z©) is -0.12 and standard deviation is 0.1. Solar 
metallicity is taken to be 0.02. 
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15 One could argue we will estimate lower limits to the scatter of the CMR 



In order to separate between the effects of age and metallicity 
scatter on the CMR scatter, we first consider the effects of for- 
mation history alone in Fig. [9] One can see a considerable color 
scatter from the ongoing accretion of recently-truncated galax- 
ies, along with some low-level bimodalities caused by kinks in 
the time evolution of colors for some metallicities (see, e.g., 
Fig. |6]l. For solar metallicity Z = 0.02 the color evolution is 
smooth, giving no 'pile-ups' in color space. The strong influ- 
ence of metallicity on the color of the red sequence can also be 
seen in this plot. 

In Fig. [10] we show our expectation for the evolution of 
the CMR intercept and scatter as a function of redshift. We 
have drawn galaxies from a log-normal distribution in metallic- 
ity with mean [Fe/H]= -0.12 and scatter 0.1 dex which is held 
constant over redshift (Fig. II lll lo i This distribution was cho- 
sen to approximately match the observed CMR intercept and 
scatter, but it is interesting to note that this value of metallic- 
ity scatter is consistent with the estimated intrinsic scatter in 
metallicity of present-day early-type galaxies of ~ 0. 1 dex from 

with this approach. Such an argument has merit, but would neglect (as we also 
have) the influence of any age/metallicity anti-correlation on the CMR scatter 
(as observed by e.g. Trager et al. 2000) 

16 Our assumption of a redshift-independent scatter is clearly an over- 
simplification; yet, in the context of our model in which each galaxy has its own 
single metallicity it is a defensible one. In reality, the metallicity may evolve 
if there is any low level residual star formation in early-type galaxies (indeed, 
if there is an age-metallicity relation in the stars in an individual galaxy, its 
light-weighted age can evolve as the galaxy ages even in the absence of star 
formation). 
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iGallazzi et al.l d2006l) . We show the mean and scatter (calcu- 
lated with the same estimation algorithm as the scatter in the 
observed relations) of the distribution as solid lines. The re- 
sulting distribution agrees well with the scatter measurements 
at the redshifts 0.05, 0.55 and 0.7. At z= 1.0, the measured 
scatter exceeds the model prediction (allowing for more scatter 
in star formation histories). The results are presented in Table 

m 

It is interesting to compare our results with a similar analysis 
of ear ly-type galaxies in clusters carried out bv lvan Dokkum & Franxl 
d200lh . They use a similar model of a constantly-growing red 
sequence to predict both the mass-to-light ratio and color evolu- 
tion (and scatter) of the early-type galaxy population, finding a 
relatively slow color evolution and a relatively constant scatter. 
Our results are in qualitative accord with their results for the 
cluster early-type galaxy population, but are representative of 
the cosmic-averaged early-type galaxy population and are con- 
strained to reproduce the observed build-up of the early-type 
galaxy population. 

5.4. Blue Spheroids 

So far, we have made no particular assumption on the mech- 
anism by which galaxies quench their star formation and join 
the red sequence. In this section, we explore a rather more spe- 
cific scenario: galaxies structurally transform into an early-type 
galaxy (e.g., through galaxy merging) before reddening on to 
the red sequence. Such a scenario predicts that there should 
be a non-negligible population of blue spheroids - i.e., galax- 
ies that are structurally early-type but have blue colors. In this 
context, we present a range of predictions for their abundance, 
requiring that blue spheroids i) have colors at least 0.25 mag 
bluer in U — V than the locus of red sequence galaxies at that 
redshift, and ii) are at least a given time t 



recognize 



from the trun- 
cation (or spheroid creation) event. This latter criterion is in- 
troduced to allow for a delay t, 



recognize 



between the event that 
truncated star formation (e.g. a galaxy merger) and the galaxy 
becoming recognizably early type; in what follows we explore 
the range < t recogn i ze < 0.4 Gyr, motivated by simulations of 
galaxy mergers. The ratio of blue spheroids to the total number 
of spheroidal galaxies (red and blue spheroids) Nbs /(Nrs+Nbs) 
is measured as a function of redshift and is shown in Fig. [T2lfor 
3 different values of trecognize- 

There are two key points to take away from Fig. [T2] First, the 
predicted blue spheroid fraction depends sensitively on one's 
choice of t recog „i ze (and therefore, also, details of the star for- 
mation history and dust content of the galaxyQ Second, 
despite this uncertainty, the range in model prediction is in 
reasonable agreement with observed blue s pheroid fractions 
(e.g. N BS /(Nrs+Nbs) ~ 0.06 at z ~ 0.6 fromlHauBleil l200f *h. 
The fraction observed by Bamfor d et alj d2008l) for the relevant 
mass range is slightly higher than in our models, but as they 
used a different method these values might not be directly com- 
parable. Recall that this toy model was constructed to explore 

17 In fact, for long t recog , n - e , the blue spheroid fraction decreases towards 
higher redshift because t recogn i ze becomes comparable to the time taken to tran- 
sition between the blue cloud and red sequence. 

18 Here we are comparing with the blue spheroid fraction of relatively mas- 
sive spheroidal gala xies, which ha ve a mass ran ge consistent w ith our red se- 
quence galaxies (see HauBler 2007, Fig. 4. 13). As HauBler 2007 (see paragraph 
4.5.1) found many of the works quoting a higher blue spheroid fraction are in- 
cluding low-mass / low-density objects tha t can never turn into a present-day 
bulge-dominated red sequence galaxy (e.g., Abraham et al. 1999; Schad e et alj 
[l999l;IIm et all200l) . 



the implications of the growth of the red sequence through the 
truncation of star formation in blue galaxies. Here we have 
shown that such a model reproduces simultaneously the evolu- 
tion of the intercept and scatter of the color-magnitude relation 
since z = 1 and the blue spheroid fraction at intermediate red- 
shift. Such an analysis lends considerable weight to the notion 
that the early-type galaxy population has grown considerably 
between redshift one and the present day through the trunca- 
tion of star formation in blue galaxies (through mergers or some 
other physical process). 

6. DISCUSSION AND CONCLUSIONS 

The evolution of the scatter of the red sequence is an im- 
portant source of information about the evolutionary history of 
the red sequence/early-type galaxy population. In this paper, 
we constructed high-accuracy color-magnitude relations at four 
different redshift ranges using accurate color and spectroscopic 
redshift information. We used a sample of over 3000 galaxies 
in the CDFS, in conjunction with a local reference sample of 
galaxies from the SDSS, to understand the redshift evolution of 
CMR scatter. We used images from the Gems survey taken 
with the ACS onboard HST to provide both high-resolution 
morphologies (to classify galaxies as early-type) and to mea- 
sure accurate colors within the half-light radius. In order to cal- 
culate accurate ^-corrections to rest-frame passbands, we used 
spectroscopic redshifts compiled from a variety of sources. At 
all redshifts, we apply a structure and color cut to isolate early- 
type red-sequence galaxies; at intermediate redshift we were 
able also to excise star-forming galaxies and AGN from the 
sample using X-ray and 24/im information. 

The resulting scatter of the color magnitude relation is ~ 0. 1 
mag in U — V color at z = 0.05, 0.55 and 0.7, and some- 
what higher at z = 1 .0. This scatter is comparable to those 
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FIG. 12. — The predicted ratio of blue spheroids to the total number of 
spheroids as a function of redshift for our toy model for the growth of the 
red sequence. The three curves show the effect of three different choices for 
trecognize < the time taken for a galaxy to become recognizably early type after 
the event that truncates star formation. Despite the dependence of Nss/(Ngs + 



Nbs) on f re 



it is remarkable that the predicted blue spheroid fraction is 



in the range of the observations (e.g. HauBler et a U2007T) . 
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fo und in some local gala xy clusters (e.g. Abell 85 and 754 
of iMdntosh etail l2005bt their Table 8) and larg er than that 
found for the Coma and Abell 496 clusters (au-y ~ 0.05 mag ; 
iBower et alJl992tlTerlevich et al.l200ll;lMcIntosn et alj20 05b). 
We note in passing that a better observational handle on the en- 
vironmental dependence of the scatter of the color-magnitude 
relation, with special attention being paid to possible cluster- 
to-cluster differences, would be highly desirable (but is clearly 
beyond the scope of this paper). 

We explored the implications of a measured CMR scatter of 
au-y ~ 0. 1 mag on our understanding of the evolution of early- 
type galaxies. The CMR scatter can be influenced by when 
star formation started, when star formation stopped, residual 
star formation and metallicity. Thus, it is not straightforward 
to turn the measured CMR scatter into robust constraints in this 
multi-dimensional space. 

We adopted instead a simple approach. The evolving number 
density of red sequence galaxies has been measured — assum- 
ing that this evolution in number density is from the quench- 
ing of star formation in galaxies with roughly constant star for- 
mation rate before quenching, we have a prediction (motivated 
by other observations) for the quenching/truncation history of 
galaxies going on to the red sequence. We have also a measure- 
ment of t he present-day meta llicity scatter of early-type galax- 
ies from Galla zzi et alj d2006l) . Thus, we have both elements in 
place where we can predict the evolution of the CMR scatter. 
Importantly, this model — in which new red sequence galaxies 
are being constantly added at the rate required by observations 
of red sequence galaxy number density — predicts the correct 
scatter in the colors of these red sequence galaxies. 

Furthermore, this model predicts approximately the correct 
number density of blue spheroids - galaxies which are struc- 
turally early-type but have blue colors - although admittedly 
with considerable model dependence. Thus, we conclude that 
these different observations - the evolution of the number den- 
sity of red sequence galaxies, the evolution of the intercept and 
scatter of the color-magnitude relation, and the blue spheroid 
fraction - paint a consistent picture in which the early-type 
galaxy population grows significantly between z = 1 and the 
present day through the quenching of star formation in blue 
cloud galaxies. 
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FIG. 13. — All objects with spectroscopic redshifts adopted in this paper are 
shown here at their positions on the sky. The different symbols showing the 
origin of the spectroscopic redshift value for each object. 



APPENDIX 

SPECTROSCOPIC REDSHIFT SAMPLE 

In order to study the evolution of the scatter of the red se- 
quence, one must use a sample of galaxies with spectroscopic 
redshifts. Toward this end, we compiled a catalog of spectro- 
scopic redshifts from a variety of sources, listed below. 

We used the data of some projects which surveyed the en- 
tire CDFS but a better coverage is reached for the Goods field 
(see also Fig. [13) . A collection of spectroscopic redshifts for 
the CDFS can be found in the 'Master Catalog of CDFS spec- 
troscopy' ( http://www.eso.org/science/goods/spectroscopy/ 
CDFS_Mastercat/). We used dat a from one paper included in 
this catalog dSzokolv et al]l200l Also the VVDS (VIMOS- 
VLT Deep Survev: lLe Fevre et alj20"04T) . lRavikumar et alJd2006l) 
and S. Koposov et al. (in preparation; VVDS high resolu- 
tion spectroscopy) surveyed the CDFS. With data limited to 
the Goods field we used the VLT/FO RS2 Sp ectroscopy Data 
Release 3.0 (IVanzella et al.ll2005l 120061 120081) a nd the GOODS 
VLT/ VIMOS Spectroscopy Data Release 1.0 dPopesso et alj 
2008). For sources with good quality flags the spectroscopic 
redshifts were taken without further testing. These are: 



VVDS - 4 (confidence level 100 per cent) and 3 (95 per 
cent) 

FORS2 - A (solid redshift determination) and B (likely 
redshift determination) 



ISzokolv et all ((2004) - 3.0 (reliable redshift 
determination with unambiguous X-ray counterpart) 
and 2.0 (reliable redshift determination) 

Ravikuma r et alj d2006l) - secure observations 



In the case of insecure or unknown data quality, the spectro- 
scopic redshifts were tested against photometric redshifts from 
the COMBO- 17 survey. In such cases we used only those 
spectroscopic redshifts which differ by less than 0. 1 from the 
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Table 3 

Spectroscopic redshifts from different catalogs 



redshift catalog 


number of objects used 


VVDS (Le Fevre et al. 2004) 


629 


Szokolv et al. (2004) 


38 


GOODS/FORS2 DR 3.0 fVanzella et al. 2008) 


464 


Ravikumar et al. (2006) 


303 


S. Koposov et al., in prep. 


510 


VIMOS DR 1.0 (Popesso et al. 2008) 


1086 



photometric ones. In total this leads to a sample of 3440 ob- 
jects (3030 useable for this project). See Figure[T3lfor the po- 
sitions on the sky of the objects used for flux measurements. 
The total numbers of objects of each catalog is listed in Ta- 
ble [3] The combined sample covers a roughly constant fraction 
of about 20% of the galaxies in the relevant magnitude range 
(20 mag < Mj < 22 mag) and then drops off towards fainter 
magnitudes with a fraction of about 10% forM/ < 24 mag. Red 
and blue galaxies show a similar behaviour. 
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